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ABSTRACT 

The abundance of compact, massive, early-type galaxies (ETGs) provides important 
constraints to galaxy formation scenarios. Thanks to the area covered, depth, excel¬ 
lent spatial resolution and seeing, the ESO Public optical Kilo Degree Survey (KiDS), 
carried out with the VLT Survey Telescope (VST), offers a unique opportunity to 
conduct a complete census of the most compact galaxies in the Universe. This paper 
presents a first census of such systems from the first 156 square degrees of KiDS. 
Our analysis relies on g-, r-, and i-band effective radii {Re), derived by fitting galaxy 
images with PSF-convolved Sersic models, high-quality photometric redshifts, Zphoti 
estimated from machine learning techniques, and stellar masses, M*, calculated from 
KiDS aperture photometry. After massiveness (M* >8 x 10^° M 0 ) and compactness 
{Re ^ 1.5 kpc in g-, r- and i-bands) criteria are applied, a visual inspection of the 
candidates plus near-infrared photometry from VIKING-DRl are used to refine our 
sample. The final catalog, to be spectroscopically confirmed, consists of 92 systems 
in the redshift range z 0.2 — 0.7. This sample, which we expect to increase by a 
factor of ten over the total survey area, represents the first attempt to select massive 
super-compact ETGs (MSCGs) in KiDS. We investigate the impact of redshift system- 
atics in the selection, finding that this seems to be a major source of contamination in 
our sample. A preliminary analysis shows that MSCGs exhibit negative internal colour 
gradients, consistent with a passive evolution of these systems. We find that the num¬ 
ber density of MSCGs is only mildly consistent with predictions from simulations at 
z > 0.2, while no such system is found at z < 0.2. 

Key words: galaxies: evolution - galaxies: general - galaxies: elliptical and lenticular, 
cD - galaxies: structure. 


1 INTRODUCTION 


The understanding of the physical processes which drive the 
galaxy mass built up and size accretion are among the most 
timely topics in galaxy evolution studies. Massive early-type 
galaxies (ETGs) are found to be much more compact in 
the past than in t h e present Universe (lDaddi_et_^;^j200^ 
iTruiillo et al.l I 2 OO 6 I : iTruiillo et al.l l2007l : Ivan der Wei et ^ 
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l2008tl . At redshifts z > 2, while the massive star-formin g 
disks have effective radii of several kpc llGenzel et al.|[200^ . 
the quenched spheroids (“red nuggets”) have small effective 
radii of about 1 kpc. Such red nuggets are thought to form 
through a chain of different processes: a) accretion-driven 
violent disc instability, b) dissipative contraction resulting 
in the formation of compact, star-forming “blue nug gets”, 
c) quenching of star formation llDekel k. Burkerl]|2014 1. Af¬ 
ter these processes occur, a gradual expansion in size of the 
red nuggets may take place, leading to the formation of the 
massive ETGs we observe in the nearby Universe. Theoret- 
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ical studies point to dry-mergers as the dominant mecha¬ 
nism for the si ze and stellar mass gro wth of very dense mas¬ 
sive galaxies ijKhochfar fc SilM l200fil i. In particular, minor 
mergers would provide a modest s tellar mass accretion, but 
a strong evolution in galaxy siz e (|yaiiJ3pldami.et__al, I I 2 OIOI: 

iHilz et SII 2 OI 3 I : iBelli et ^l2014l : iTbrtora et al.ll2014l l. Merg¬ 

ers are believed to be common for very massive systems at 
high redshifts, with major merger rates (mergers per galaxy 
per Gyr) in the range 0.3 — 1 Gyr~^ at z ^ 2 a nd smaller 
than 0.2Gyr“^ at z < 0.5 llHopikins et al.|[2010ll . An alter¬ 
native scenario explains the size evolution as the result of 
(e.g.) quasar feedback, rather than merging, making galax- 
ies to puff up after a lost of large amounts of (cold) gas 
dFan et al.ll2008l. l201(]h . 


Over cosmic time, one may expect that high-z 
compact galaxies evolve into present-day, massive, big 
galaxies. However, a fraction of these objects might 
survive intact till the present epoch, resulting in com¬ 
pact, relic systems in the nearby Universe characterized 
by old stellar populations. Recently there have been 
some efforts to sear ch for massi v e compact galax¬ 
ies a t low redshifts _( [Truiillo et ajj_j2009t__lTavlor et al.l 
2010|;__^alentmuzzie£_alJ_[20]U|;__lShih^_&Stockton|_[201ll 


^ujiUoe£_alJj2012jPoeEiantiet_alJ 20WaHPoggianti_^^^ 


2013bl : iDamianov et all l20Rj . |2014| . l2015al : IS^lder etld 

2 OI 5 II and investigate further their dynamical and stellar 
population properties , as well as the role of environment on 
their pr operties (e.g. Ferre- Mat^ et ahl 201^ Laskeret^^y 
2 OI 3 I: iFerre-Mateu et al.l l2015l: Thuiillpet_^y__l20lJ; 


, Damianov et al.l 2015bl : Stringer et al.l 
l2015ali. 
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In some theoretical models, that include the effect of 
galaxy mergers, the fraction of massive objects that would 
survive without undergoing any significant transformation 
since z ~ 2 to the present could reach a fraction of a bout 
1 — 10% (iHopkins et al.l [20091 : lOuilis fc Truiillol l2013l l. At 
“low” redshifts (z < 0.2), theoretical models predict a den¬ 
sity of relic remnants in the range 10“^ — 10“®, which 
means that, in large surveys like the Sloan Digital Sky 
Survey (SDSS), we might expect to find a few candidates 
in this redshift range. However, recent observational works 
have shown the paucity of old, super-compact (Re < 1.5 
kpc), massive (M ^ >10^^ M(:^) galaxi e s in the loca l Uni¬ 
verse (z < 0.2; iThdillo et ahl^OOgl : iTavlor et al.l I 2 OIOI : 
iFerre-Mat^ et al.ll2012l l. Indeed. NGG1277, a nearby lentic¬ 
ular galaxy in the Perseus cluster, is actually the only well- 
charact erized, old system at z ~ 0, th at might be a true relic 
galax y (iTruiillo et al.l [2014I : see also iMartin-Navarro et al.l 
2 OI 5 II. Other candi dates have been recently detected by 
Saulder et al.l ll2015l l. although only a few of them fulfill the 
above (restrictive) size and mass criteria (Re < 1.5 kpc and 
M* ^10^^ Mq), none of them being at z < 0.05. Larger num¬ 
bers of old compact systems have been found at lower masses 
(< 10 ^^ Mp)), when relaxing the compactness selectio n cri¬ 
teria dValentinuzzi et ^l2010l : IPoggianti et al.]l2013al l. The 
(almost) lack of nearby super-compact, relic systems may 
represent a challenge for the current paradigm of galaxy for¬ 
mation; in particular, we have to understand whether this 
lack is due to some observational bias, e.g. the limited spa¬ 
tial resolution of photometric data at z ~ 0; a failure of 
theoretical predictions; and/or an environmental effect, for 


which, as suggested by NGC 1277, relic galaxies might be 
more frequent in high-density cluster regions. 


In the intermediate redshift range (0.2 z ^ 0.7), 
compacts have been recently investigated in detail by 
IDamianov et al.l ll2014fl . who selected ~ 200 massive com¬ 
pacts from a sample of stellar-like objects within the 
6373.2 sq. deg. of the BOSS survey; 20% of these galaxies 
are dominated by old stellar populations, which make them 
reliable candidates to be the product of the unperturbed 
evolution of compact high-z systems. However, 93% of these 
galaxies do not have measured Re which hampers the selec- 
tion of such systems as tr uly compact objects. More recently, 
IDamianov et al.l ll2015al L have analyzed F814W HST im¬ 
ages for the COSMOS field, providing robust size measure¬ 
ments for a sample of 1599 compact systems in the range 
0.2 < z < 0.8. Other studies have performed detailed anal¬ 
ysis of stellar populations and morpholog y of small sam¬ 
ples of compact galax ies at these redshifts I Hsu et aJll2014l : 
IStockton et al.l l2014h . The population of dense passively 
evolving galaxies in this intermediate redshift range possibly 
represents a link between compact systems, dominating the 
massive quiescent galaxy population at high z, and their de¬ 
scendants in the nearby Universe. Indeed, large samples of 
compacts, with high-quality photometry (to derive reliable 
structural parameters), and spectroscopic data, are actually 
necessary to better understand the formation and evolution 
of these systems. 


The Kilo Degree Survey (KiDS; Ide Jong et al. I I 2 OI 5 I) is 
one of the ESO publi c surveys carried out with the VLT Sur¬ 
vey Telescope fVST: ICapaccioli fc SchiDanill201lh equipped 
with the one square degree field of view an d high angular 
resolution (Q.2"/pix el) OmegaCAM camera llKuiiken et al.l 
I 2 OO 4 I : iKuiikenll^Ol il l . KiDS is mainly designed for weak leas¬ 
ing studies, providing deep imaging in four optical bands 
(ugri), over a 1500 square degree of the sky with excellent 
seeing (e.g. 0.65” median FWHM in r-band). The high im¬ 
age quality and deep photometry are ideal to investigate 
massive compact systems. 


_ A ccordi ng to predictions from simulations iIGuo et al.l 

l2011al. l20ldL we can expect to find ~ 0.3 — 3.5 relic per 
square degree, at redshift z < 0.5. This prediction does criti¬ 
cally depend on the physical processes shaping size and mass 
evolution of galaxies, such as the relative importance of ma¬ 
jor and minor galaxy merging. 


Several compa ct galaxy definition s have been adopted 
in the literatur e llTruiillo et al. 200^ Taylor et al.l I 2 OI 0 I : 
IPoggianti et al.l l2013al : IDamianov et al.l l2015al L In the 
present work, we present the properties of a sample of dense 
massive galaxy candidates in KiDS, defining as massive 
super-compact galaxies (MSCGs) those early -type systems 
with M* > 8 X 10^° Mq and Re < 1.5 kpc (ITruiillo et al.l 
l2009tl . These selection criteria are rather conservative, pro¬ 
viding the ideal benchmark for galaxy evolution theories. 
The paper is organized as follows. In Sec. [2] we present the 
KiDS data sample and the selection of our photometrically 
selected compact galaxies. The main results, and in partic¬ 
ular the evolution of number density as a function of red¬ 
shift, are presented in Sec. [S] A discussion of the results 
and future prospects are outlined in Sec. U) We adopt a cos¬ 
mological model with (Dm, DA,h) = (0.3,0.7,0. 70), where 
h = 77o/100 kms“^ Mpc“^ llKomatsu et al.ll2011^ . 
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2 SAMPLE SELECTION 

The galcixy sample presented in this work is based on the 
data include d in the first and sec ond data releases of KiDS 
presented in ide Jong et al.l (l2015ll . which we address the in¬ 
terested reader for details. The total dataset includes 156 
KiDS pointings (133 from the KiDS data release 2), in which 
we have identified about 22 million sources including ~7 
millions which have been classified as high quality extended 
sources (mostly galaxy-like, see below). A full description 
of the galaxy sample is given in Napolitano et al. (2016, 
in prep.). In the following section, we summarise the main 
steps for the galaxy selection procedure and the determina¬ 
tion of galaxy physical quantities as structural parameters, 
photometric redshifts and stellar masses. 


2.1 KiDS high signal-to-noise galaxy sample 

We start from the KiDS multi-band source catalogs, 
where the photometry h as been obtained with S-Extractor 
dSertin fc ArnoumllOOd ) in dual image mode, using as ref¬ 
erence the positions of the sources detected in the r-band 
images. While magnitudes are measured for all of the filters, 
the star/galaxy separation, positional and shape parame¬ 
ters are based on the r-band data. The choice of r-band is 
motivated by the fact that it typically has the best image 
quality and thus provides the most reliable source positions 
and shapes. Critical areas as saturated pixels, star spikes, re¬ 
flection halos, satellite tracks, etc. have been masked using 
both a dedicated automatic procedure and visual inspection. 
The total area after removing the borders of the individual 
KiDS pointings is 163.9 square degrees, while the unmasked 
effective area adopted is of 105.4 square degrees. 

Star/galaxy separation is based on the distribution of 
the S-Extractor parameters CLASS_STAR a nd S/N (signal-to- 
noise rati o) of a number of su re stars Isee lLa Barbera et al.l 
I2OO8I and ide Jong et ^l2015l for further details). We have 
further retained those sources which were marked as be¬ 
ing out of critical area from our masking procedure. From 
the original catalog, the star/galaxy separation leaves ~ 11 
million galaxies, of which ~ 7 millions have high quality 
photometry being non-deblended sources located out of the 
masked area. 

To perform accurate surface photometry and determine 
reliable structural paramete rs, the highest-quality source s 
have been further selected dLl Barbera et al.l l2008l . I2OI0II . 
Thus, we have finally gathered those systems with the high¬ 
est S/N in the r-band ima ges, S /N_^ = l/MAGERR_AUTO_r> 
50 dLa Barbera et ^ ] |2008l . [2010l: Roy et al. 2016, in prepa¬ 
ration). This sample consists of ~ 380 000 galaxies. 

Relevant properties for each galaxy are derived as de¬ 
scribed here below: 

• Photometry. As a standard KiDS catalogs’ parameters, 
we have derived S-Extractor aperture photometry in the 
four bands (ugri) within several radii. For our analysis we 
have adopted aperture magnitudes MAGAP_2, MAGAP_4 and 
MAGAP_6, measured within circular apertures of 2", 4” and 
6 ” of diameter, respectively, while a first probe of the total 
magnitude is provided by the Kron-like MAG_AUT0. 

• Structural parameters. Surface photometry is performed 
using the 2DPHOT environment. 2DPHOT is an automatic 


software designed to obtain both integrated and surface pho¬ 
tometry of galaxies in wide-held images. The software hrst 
produces a local PSF model from a series of identihed sure 
stars. This is done, for each galaxy, by htting the four clos¬ 
est stars to that galaxy with a sum of three two-dimensional 
Moffat functions. Then galaxy snapshots are htted with 
PSF-convolved Sersic models having elliptical isopho tes plus 
a local background value fsee lLa Barbera et ahlfioo^ for fur¬ 
ther details). The ht provides the following parameters for 
the four wavebands: surface brightness ^e, circularized effec¬ 
tive radius, Re, Sersic index, n, total magnitude, ms, axis 
ratio, q, and position angle. As it is common use in the liter¬ 
ature, in the paper we use the circularized effective radius. 
Re, dehned as Re = ^/qRe,maj, where Re.maj is the major- 
axis effective radius. 

• Photometric redshifts. Redshifts are determined with 
the Multi Layer Perce ptron with Quasi Newton Algorithm 
fMLPQNAl method (IBr escia et alJl20lA l2014l l. and fully 
presented in ICavuoti et al.l ll2015ll . which we refer for all de¬ 
tails. Both apertures of 4" and 6" of diameter are used. In 
machine learning supervised methods, a knowledge sample 
is needed to train the neural network performing the map¬ 
ping between magnitudes and redshift. The knowledge base 
consisted on spectroscopic redshift fr om the Sloan D igital 
Sky Survey data release 9 fSDSS-DR9: lAhn et al.l[2012|j and 
Galaxy And Mass Assembly data release 2 (GAMA-DR2; 
[Driver et al.l I 2 OI 1 I I which together provide redshifts up to 
z ^ 0.8. This knowledge base includes ~ 60000 objects, 60% 
of which are used as training sample, and the remaining ones 
are used for the blind test set. The redshifts in the blind test 
sample resemble the spectroscopic redshifts quite well. The 
scatter in the quantity Az = (zapee — 2phot)/(l + 2spec) is 
~ 0.03. After these experiments, we have finally produced 
the final catalogue of redshifts for our sample. The cut op¬ 
erated in the fitting magnitudes to resemble the luminosity 
ranges in the knowledge base will impact the completeness 
in the faint regime. But our high-5/A sample is not affected 
by this source of incompleteness. 

• Stellar masses. We have used the software Le Phare 
llArnouts et al.l 1 19991 : lllbert et al.l l2006lj , which performs a 
simple fitting method between the stellar population syn¬ 
thesis (S PS) theoretical mod e ls and data. S ingle burs t mod - 
els from iBruzual fc Gharlotl (l2003l l and a IChabrieil ll200l[ l 
IMF are used. Models are redshifted using the photomet¬ 
ric redshifts. We adopt the observed ugri magnitudes (and 
related 1 a uncertainties) within a 6” aperture of diame¬ 
ter, wh ich are corrected for Galact ic extinction using the 
map in lSchlaflv fc Finkbeineil 1I2OIIII . Total magnitudes de¬ 
rived from the Sersic fitting, ms, are used to correct the 
outcomes of Le Phare for missing flux. The single burst as¬ 
sumption is suitable to describe the old stell ar populations 
in the compact galaxies w e are interested in (iThomas et al.l 
I2OO5I : iTortora et al.ll200^ . The estimated photometric ages 
are used to check if galaxies are compatible with being relic 
remnants of systems formed at 2 ~ 2. The degeneracy be¬ 
tween age and metallicity in the stellar population analysis 
can be solved only with forthcoming spectroscopic follow¬ 
ups. 

• "Galaxy classification". Using Le Phare, we have also 
fitted the observed magnitudes MAGAP_6 with a s et of 66 
empirical spectral templates used in lllbert et al.l 1I2OO6II . 
The set is based on the four basic templates (Ell, Sbc, 
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masses ^5 x 10^° up to redshift ^ = 0.5 (see see Napolitano 
et al. 2016, in prep., for further details). 



J-K 


Figure X. J — K vs. g — J diagram. An aperture of 2" of di¬ 
ameter is adopted here. Blue and cyan symbols are for high- 
confidence stars with any error and 5J = 5K < 0.05, respectively. 
Red points and open boxes are for compact candidates with any 
error and SJ = SK < 0.05, survived after the criteria on the 
mass, size, 2D fit quality and visual inspection. Vega J and K 
magnitudes are converted to AB as Jab = -^Vega + 0.930 and 
Kab = -^Vega + 1.834. We highlight the regions which are popu¬ 
lated by stars (blue), galaxies (yellow) and QSOs (pink). We have 
considered as sure galaxies those objects with colours J — K > 0.2 
and g — j >2 (yellow shaded region). 


Scd, Irr) des cribed inIColeman et al.l (|l980h . and star burst 
models from Kinney_eL_^ ( 19961) . GISSEL synthetic models 
(|Bruzual &; Chariot! |2Q0^ are used to linearly extrapolate 
this set of templates into ultraviolet and near-infrared. The 
final set of 66 templates (22 for ellipticals, 17 for She, 12 for 
Scd, 11 for Im, and 4 for starburst) is obtained by linearly 
interpolating the original templates, in order to improve the 
sampling of the colour space. We have selected the ETGs 
by choosing the galaxies which are best fitted by one of the 
elliptical templates (see Napolitano et al. 2016, in prepara¬ 
tion). 

• Near Infrared photometry from VIKING-DRl. As a 
complementary dataset for our selection of compact galax¬ 
ies, we have used the data from the first release (DRl) of 
the ESO VISTA Kilo Degree Infrared Galaxy (VIKING) 
survey. The VIKING survey is the KiDS twin survey, and 
provides the near-IR coverage of the same sky region in the 
5 wavebands Z, Y, J, H and K s. VIKING-DRl consists of 
108 observed tiles (lEdee et al.ll20lj ). We have found that it 
overlaps with ~ 58% of our 156 sq. deg. in the KiDS survey. 
For uniformity, we have used photometry within the same 
aperture as the optical data: APERMAG_3 for VIKING and 
MAGAP_2 for KiDS, which correspond to an aperture of 2" of 
diameter. In particular, we have concentrated on J and Ks 
passbands, which we do not use for SPS fitti ng, but only to 
improve the star-galaxy s eparation criterion llMaddox et al.l 
l2008l: iMuzzin et al.|[20ll) as will be discussed in details in 
the next section. 


The sample of high-S/V galaxies is complete down to a 
magnitude of MAG_AUT0_r ~ 21, which correspond to stellar 


2.2 Selection of compact galaxies 

MSCGs have been selected using the following criteria: 

(i) Massiveness. The most massive galax ies with M* > 
8 X 10^'^ Mq are taken ll Trujillo et al. 1 120091 ). reducing the 
original sample of ~ 380000 galaxies to ~ 30600 massive 
galaxies. 

(ii) Compactness. We select the densest galaxies by fol¬ 
lowing recent literature l|Truiillo et al.lfioODll . We get galax¬ 
ies with median circularized radius, Re, among the g-, r- and 
i-bands, less than 1.5 kpc. About 1300 compact candidates 
remain after this selection. 

(iii) Best-fitted effective radii. The best htted structural 
parameters are considered, taking those systems with a re- 
duce d Y?n from 2DPH0T sm aller than 1.5 in g, r and i fil¬ 
ters llLa Barbera et ^11201011 . To avoid any accidental wrong 
fits we have also removed galaxies with unreasonable best- 
fitted parameters, applying a minimum value for the size 
{Re = 0.05 arcsec), the Sersic index (n > 0.5) and the axial 
ratio {q = 0.1) in all the bands. The minimum value in the 
Sersic index has also allowed us to possibly remove misclas- 
sified stars, which are expected to be fitted by a box-like 
profile (mimicked by a Sersic profile with n —>■ 0). These last 
criteria on the quality of the structural parameters reduce 
further the sample to the 106 highly reliable candidates. 

(iv) Eye-ball check. We have made an eyeball inspection 
of the images and residuals from Sersic Htting of the can¬ 
didates, with the aim of removing problematic objects or 
possible misclassihed stars. To reduce subjectivity, three 
of the authors have independently checked the images and 
graded them according to the following scheme: grade-2 are 
sure galaxies, grade-1 are uncertain galaxies lacking a well- 
defined elliptical shape, and grade-0 are misclassihed objects 
(either stars or corrupted Hts). The mean of the three classi- 
hcations has been adopted as the hnal grade. In this way, we 
have retained those galaxies with a grade larger than 1, to 
include only systems graded by at least one of the observers 
as a sure galaxy. The candidates are further reduced to 96, 
after removing objects with signiheant contamination from 
neighbours. 

(v) Optical-hNIR star-galaxy separation. We have 
adopted a mor phological criterion t o perf o rm the star-galaxy 
classi heation llBertin &; Arnouta 1 19961 : iLa Barbera et ahl 
I 2 OO 8 II and used visual inspection as ultimate check of the 
galaxy classiheation. However, based on optical data only 
a star can be still misclassihed as a galaxy on the basis 
of its morphology, and this issue can be highly dramatic 
for very compact objects (generally with size comparable 
or smaller than the seeing). In absence of spectroscopic 
information, optical-fNIR colour-colour diagrams can 
provide a stro ng constraint on the galaxy nature of the 
candidates fe.g. lMaddox et ^l2008l : iMuzzin et al.ll2013l') . In 
particular, g, J and Ks magnitudes within 2" of diameter 
are adopted for this purpose on those helds with coverage 
by the two surveys. The Vega VISTA magnitudes are 
converted to AB using the conversion formulae from Le 
Phare, Jab = Jvega + 0.930 and Kab = Kvega + 1.834. 
Stars and galaxies with the best J and Ks photometry are 
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Figure 2. 2D fit output for 4 example candidates from 2DPHOT procedure. For each galaxy, the top panels show the galaxy image 
(left) and the residual after the fit (right), while the six bottom panels provide residuals after subtraction, plotted as a function of the 
distance to the galaxy center, with each panel corresponding to a different bin of the polar angle. Residuals are normalized with the 
noise expected from the model in each pixel. 


also considered (JJ, SK < 0.05). The results are shown in 
Fig.m Stars have blue J — K colours (i.e., J — K ^ 0.2, 
see light blue shaded region in Fig. [T]). We have also found 
four objects with J — K < 0.2 and, indeed, are erroneously 
classified as galaxies. We take as compact candidates those 
systems with J — K > 0.2 and g — J > 2 (see light-yellow 
shaded region in Fig. [T]). For those galaxies with available 
VIKING photometry, this selection allows our set of can¬ 
didates from the previous classification steps to be refined. 
Out of 46 galaxies with VIKING photometry, we remain 
with 42 high-confidence candidates. 

In Fig. [T]we have also highlighted the locus populated by 
point sources with red J — K { ^0.2), but blue g — J (< 2) 
colours, classi fied as stars by S-Ex tractor, which are presum¬ 
ably quasars (iMaddox et al.l[2008l . see pink region in Fig. [T]). 
However, the analysis of this class of objects is beyond the 
scope of the paper. 

To perform an homogeneous comparison of the sample 
of our compact candidates with a sample of "normal" galax¬ 
ies, using the same original catalog described in Sec. [H we 
have created two control samples (C-Samples) of galaxies 
with similar stellar masses (M* > 8 x 10^° Mq), the same 
lower limits on Sersic index [n > 0.5), axis ratio {q > 0.1) 
and Rs (J?e > 0.05 arcsec) as those for the compact systems, 
but without applying any compactness criterion. The first 
FULL C-Sample consists of all galaxies, with no restriction on 
the galaxy type, while the ETGs C-Sample consists of ETGs, 
classified as described in Sec. [21 


3 RESULTS 

3.1 The final sample 

We start by summarizing the sample we are left with for our 
analysis: after the first three criteria we have selected 106 
candidates, which are reduced to 96 after the eye-ball check 
in (iv). The matching with VIKING included 46 candidates 
and, after applying the criterion (v), we are left with a sam¬ 
ple of 42 out 46 galaxy candidates (i.e. 4/46 objects are likely 
stars according to their optical-NIR colours, corresponding 
to a ~ 10% contamination). Updating the total number, we 
are left with 92 candidates, of which ~ 4 further sources 
from the region that does not overlap with VIKING might 


be stars. This corresponds to a density of ~ 0.9 compact 
galaxies per square degree, cumulatively for z < 0.7, while 
we do not find objects at z ^ 0.2. In contrast. [Trujillo et al.l 
ll2009li find 29 secure MSCGs at z < 0.2 within the SDSS 
area, using our same criteria. However, only one of these is 
at very low redshift (z < 0 . 1 ), and none of th em features old 
stellar populations llFerre-Mateu et ^l2012|j . and thus they 
cannot actually be classified as relics. The absence of com¬ 
pact systems at z < 0.2 in our datasample might be related 
to environmental issues, since in the current KiDS area we 
are missing nearby clusters of galaxies , where a larger frac - 
tion of compact galaxies is predictedjStringe£^^^l]j 201 ^, 
and seems to be obser ved indeed llPoggianti et al.l l2013al : 
IValentinuzzi et ^ |2010|) . A simi l ar env ironment effect would 
be present in the Trujillo et al. I ll2009ll sample. Examples of 
the 2D fitting results for four candidates are presented in 
Fig. (21 where postage stamps of the galaxy images and the 
residuals are shown. 

Fig. [3]shows the distributions of redshift, axis ratio, Ser¬ 
sic index, effective radius, stellar mass and surface mass den¬ 
sity for our sample of compact candidates. The median red¬ 
shift of the sample is ~ 0.44, with an rms of ~ 0.18, larger 
than the median redshifts of the FULL and ETGs C-Samples 
(~ 0.36 and ~ 0.34, respectively). Gompact candidates have 
a median r-band effective radius of ~ 1.2 kpc (rms=0.31 
kpc), a Sersic index of ~ 4.3 (rms=2.3), an axis ratio of 
0.40 (rms=0.14) and median stellar mass of ~ 10.99 dex 
(rms=0.09). On the contrary, for the FULL C-Sample the 
median size, Sersic index, axis ratio and M* are ~ 7.7 kpc, 
4.7, 0.73 and ~ 11.07 dex, respectively. Finally, if the ETGs 
C-Sample is considered, the median size, Sersic index, axis 
ratio and M* are ~ 8.4 kpc, 5.5, 0.75 and ~ 11.10 dex, re¬ 
spectively. Thus, smaller sizes translate into shallower light 
profiles (smaller n) when compacts are compared with nor¬ 
mal gal axies (consistent ly with the galaxy merging frame¬ 
work in iHilz et al.ll20i^ . Gompact candidates also present 
lower axis ratios than normal galaxies. These elongated 


high-z (Ivan der Wei et al. 20 111: 

Buitrago et al.ll2013ll. and 

at lower redshifts iTruiillo et al.l 

20ll 2014ll. 


In our compact sample, the average surface mass den¬ 
sity within 1 i?e is ~ 1.2 x 10^° Mq/ kpc^ (rms of ~ 5.6 x 
lO®M 0 /kpc^), about 2 orders of magnitude larger than 
the average surface mass density of the C-Samples. This 
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Figure 3. Distribution of some galaxy parameters. From left to right we show: 1) photometric redshifts, 2) r-band axis ratio, 3) r-band 
Sersic index, 4) r-band effective radius, 5) Chabrier IMF-based stellar mass, 6 ) r-band effective surface mass density. 


median density is a bit larger than the range of values 
of 7 — 15 X 10^ Mct) /kpc^ found in 4 compact galaxies at 
^ < 0.2 by iTruiillo et all (|2Q12l'l . but the values are con- 
sistent within t he sca tter. Taking the only 3 gala:xies from 
ISzomoru et al. I ll2012ll . which satisfy the same criteria we 
are adopting in this paper, we find a median density of 
~ 1.1 X lO^°M 0 /kpc^, which is fully consistent with our 
result. 

The r-band Re as a function of M* is shown in Fig. ID 
The compact candidates are plotted with a random subset 
of ~ 2500 galaxies from the FULL C-Sample. As the com¬ 
pact candidates are selected to have small sizes, they lie 
in a region of the size-mass diagram where very few ob¬ 
jects are detected, which provides a further evidence in fa- 
vor of the rarity of these compact objects at z ^ 0.7 (Fig. IH 
ITruiillo et al.l 2009h. In p e ;. I^ we also overplot the 2 < 0.2 
ga laxies in TVujilloe£_^ (120120 and the three 2 ~ 2 galaxies 
in ISzomoru et al.l ( 2012 l l which fullfil our selection criteria. 
The evolution with redshift of the sizes, magnitudes and 
colours are shown in Fig. [S] The ReS look quite constant 
as a function of redshift for the compact systems. At all 
the redshifts, almost all the compacts have faint MAG_AUT0 
when compared with the control sample population (middle 
panel in Fig. [S|. We find that 75 candidates (~ 82% of the 
whole sample) lie below the r-band magnitude completeness 
limit of ~ 21 , implying that our sample is complete up to 
redshift 2 ~ 0.5. Finally, most of the galaxies populate the 
red-sequence, and are the best candidates to be remnants of 
high-z red nuggets (right panel in Fig. [5]). 

58 out of 92 galaxies have old ages consistent with a 
formation redshift 2t ^ 2 , and so could b e the remnants 
of the compact galaxies observe d at 2 > 1 l|Gargiulo et al.l 
l201ll. I 2 OI 2 I: ISzonioru et ^l2012ll . 82 out of 92 galaxies (i.e. 
~ 89% of the whole sample) are best fitted by an ellipti¬ 
cal template, and are classified as ETGs (see Sec. [2]). The 
rest of the galaxies have redshifts 2 > 0.5 and have colours 
which are best-fitted by Sbc models. However, particularly 
at these redshifts, one should take this colour-based classifi¬ 
cation with caution, as a spectroscopic follow-up is actually 
needed to perform an accurate stellar population analysis. 


3.2 Systematics from wrong redshifts 

We have finally cross-matched our sample of candidates with 
SDSS and GAMA spectroscopy catalogs, finding spectra for 
9 of the selected candidates. 

One of the main systematics in our selection of com¬ 
pact galaxies is induced by wrong redshift determinations, 
which can affect both the (linear) effective radii and stel¬ 
lar masses, moving the compact out of our selection crite¬ 
ria (e.g., see Figs. |4] and [Sj. 7 out 9 of these systems are 


included in the SDSS-fGAMA datasample used for photo¬ 
metric redshift determinations (see Sec. [2|: 5 in the training 
sample and 2 in the test sample. Although the photometric 
redshifts are shown to ap proximate quite we ll the spectro¬ 
scopic ones (see Sec. [2]and [C^avuoti et al.l201^ 1. we note that 
almost all of the photometric redshifts are underestimated 
with respect to the spectroscopic value. The median differ¬ 
ence A 2 = 2 phot — 2 :spec is —0.07, which increases to —0.1 if 
the galaxies not in the training sample are considered. This 
systematics can be probably related to the under-sampling 
of this galaxy population, just 5 galaxies in the training sam¬ 
ple, which can fail the optimum training of the network, and 
the quality of photometric redshift outcomes is degraded. 

As first test, we have re-computed the sizes and stellar 
masses of these 9 objects using the 2 spec. Only 4 galaxies 
survive to the size and mass selection criteria, i.e. ~ 44% (1 
not in the training sample). 

Despite the small sample available for this test, we can¬ 
not exclude the presence of a systematics of the photometric 
redshift determination for compact sources which we might 
quantify of the order of A 2 = —0.1. Qualitatively, if the 
spectroscopic redshift is larger than photometric one, then 
Re in physical scale gets larger. The effect on M* is not as 
simple as the one on Re- As an exercise, we have investi¬ 
gated the impact of a redshift error of A 2 = 0.1 on the 83 
galaxies in the sample without spectroscopic redshift. Using 
re-calculated values for Re and stellar mass, we find that 26 
of the sources (31%) survive the cuts. 

Therefore, a spectroscopic follow-up will be needed to 
finally confirm the nature of these galaxies, allowing us to 
classify these candidates as compact or very-compact galax¬ 
ies. Increasing the knowledge base and the population of 
compact systems with measured spectroscopic redshifts will 
improve the quality of the photometric redshift estimates. 


3.3 Colour gradients 

A big improvement of our analysis with respect to previous 
works on compact ga laxies at similar redshift ranges (e.g. 
IPamianov et al.ll2014[ l is the high-S/A^ photometry which 
allows us to derive robust structural parameters and ob¬ 
tain the first determinatio n of colour pro f iles fo r these sys¬ 
tems. Therefore, following iTortora et al.l (I 2 OI 0 I I. we define 
the colour gradient as the angular coefficient of the rela- 

c / Y" _ Yz-N 

tion X — Y vs. log 77, Vx-y ~ —7 - measured in 

0 log R 

mag/dex (omitted in the following unless needed for clar¬ 
ity). The fit of each color profile is performed in the range 
Ri = 0.1i?e ^ R ^ R2 = Re, where the effective ra¬ 
dius is measured in the r-band. By definition, a positive 
CG, Vx-Y > 0, means that a galaxy is redder as R in- 
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Figure 5. We show the median of the Re in the g-, r- and i- bands (left), r-band MAG_AUT0 (middle) and g — r (right) vs. redshift. The 
gray region in the left panel sets the region within which we have searched the compacts. Gray line in the middle panel sets the 90% 
completeness limit of the high-5/Af sample. The g — r colour is calculated within &' of diameter. See further details about the symbols 
in Fig. |4] 



totally related to the larger uncertainties on the structural 
parameters in our small objects. 


These results agree with previously reported ranges of 
color gradients of pas sively evolving massive galaxies at 


low- 

or intermediate-z (iTamura et al.ll200C 

: Tamura & Ohta 

200C 

; IWu et al.ll2005l: Tortora et al.l 2010l: 

La Barbera et al. 

2012 

1 or at high-z IIGuo et al.ll2011b: Gargiulo et al.ll201ll. 

2012 

1 and with simulations (e.g.. Tortora et al. 20131). 


Finally, we find that 30 compact candidates have all 
negative colour gradients, 9 have positive gradients, while 
the rest of the sample have at least one of the three gradients 
with a different sign with respect to the others. This wide 
range of behaviours demonstrate that these objects can be 
form e d in quite different in itial conditions (|Gargiulo et al.l 
I 2 OI 2 I : IPamianov et al.ll20l4) . 


Figure 4. We show the median of the Re in the g-, r- and 
i- bands vs. stellar mass. Compact candidates (plotted as red 
points) are compared with a selection of galaxies (~ 2500 ran¬ 
domly extracted) from the FULL C-Sample (small black points). 
Gray shaded region highlight the region where candidates are 
selected (see text for details about the adopted selection crite¬ 
ria) . We also plot the change if the spectroscopic redshift is used 
(dar k-green lines). Finally , the four compact galaxies with z < 0.2 
jromJlTujiUoet ajj_ ^^ and the three compacts at 2 ~ 2 from 
ISzomoruet‘an ' (l2012l l. ttiat fulfil our selection criteria, are shown 
as cyan stars and blue squares, respectively. 


creases, while it is bluer outward for a negative gradient. 
PSF-convolved Sersic profiles in g-, r- and i-band are used 
and, in particular, g — r,r — i and g — i colours are discussed. 
We omit detailed analysis in terms of redshift and stellar 
mass, but we pinpoint what is the range of colour gradients 
in our gala:xy sample, comparing these results with the ETGs 
C-Sample and some previous literature. 

On average, compact population has Vg_r = —0.21 
(rms 0.52), Vr-i = —0.07 (rms 0.59) and Vg_i = —0.30 (rms 
0.73) which are substantially consistent with the gradients 
for the control population of ETGs which are Vg_r = —0.17 
(rms 0.33), Vr-i = —0.05 (rms 0.23) and Vg_i = —0.21 
(rms 0.40). Hence, compact candidates look quite similar to 
normal ETGs, except for the scatter, which is partially or 


3.4 Abundance vs. redshift 


The number density of compact massive galaxies as 
a function of redshift is an important constraint on 
models of galaxy assembly. In recent years there have 
been different efforts to produce a census of such sys¬ 
tems in diff e rent re dshift bins (e . g. Tnrijllp_et_al, 200^ 
Tavlor et al.l lioiol: iGassata et al.l 20131: Poggianti et al.l 


2013bl lal: IPamianov et al.ll2014l . l2015al . Saulder^ al.l 2015h . 

Sample size is a crucial aspect to increase the constrain¬ 
ing power. If the compact galaxies found in the present 
work are a representative subsample of the whole popula¬ 
tion of compacts over the whole area of 1500 sq. deg. that 
will be mapped by the KiDS survey, we expect to increase 
the present sample by a factor of ten in the next few years. 

For what concerns our current sample, we have binned 
galaxies with respect to redshift and normalized to the co¬ 
moving volume corresponding to the observed KiDS sky 
are^B- The errors on number counts take into account fluc¬ 
tuations due to Poisson noise, as well as those due to 


^ This is obtained by multiplying the number of candidates by 
/area = TskyMsurvey, where Tgky (= 41253 sq. deg.) is the full 
sky area and ^survey = 105.4 is the effective KiDS-DR2 area. 
Then, the density is derived by dividing for the comoving volume 
corresponding to each redshift bin. 
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- KiDS 

KiDS-z_syst 

---- BOSS-Md,„ >5 X 10 '“Mo 
-- BOSS-Md|,„>8x10'“Mo 

- BOSS-M*„>10xip'“Ms 

- - BOSS - Md,„ > 6 X 1 0 '“Mo (zi > 2) 

- BOSS - Md,„ > 8 X 1 0 '“Mo (zi > 2) 

A BOSS - Re < 1.5 kpc 
COSMOS-Re < 1.5 kpc 

• Trujillo+09 

• Taylor+10 

• Poggianti+13 - Re < 1.5 kpc 
A Trijjillo+14 

■ □ Saulder+15- Re < 1.5 kpc, M. > 8 x 10'°Mo 

. Guo+11 - AM < 30% 

Guo+11 - AM < 10% 

- Guo+13 - AM < 30% 

Guo+13 - AM < 10% 

- + - Wellons+15 (lllustris) 


Figure 6. Number density of compact galaxies vs. redshift. Error bars denote Icr uncertainties (see the text). Vertical dashed black 
lines mark the completeness limit of the sample fsee lde Jong et al]|2015h . Left. Number densities for all galaxies, independent of their 
photometric formation redshift, are plotted with open squares and solid black lines. The solid gray line with open squares takes into 
account possible systematics in the redshift deter mination (see Se c. Short-dashed, long-dashed and solid violet lines are number 

densities of stellar-like objects from BOSS-DRIO ijPamianov et ani2014h with M^yn > 5 x 10^*^, 8 x 10^*^, and 10^^ LIq, respectively. 
The violet shaded region is the Icr error for the case with M^yn > 10^ ^ The viole t trian gle with bars also shows the abundance of 

galaxies at z 0.25, with Re < 1.5 kpc and M^yn > 8 x 10^*^ M(7), from lDamj anov et ^ |[2o3)* The yellow line with lighter yellow region 
plot abundances for compacts in the COSMOS area llDamianov et al1]20^3j . selected with the same criteria as in the present work 
(M* > 8 X 1 0^^ Mq) and Re < 1.5 kpc; I. Damjanov, p rivate communication ) . Red, cyan and green points are the results for compact 
galaxies from PlVuiillo et alj ||2009|). [Taylor et al.| l|2 01oh and IPoggianti et alj ll2013al) . respectively. Orange boxes show the abundances 
for compacts in SDSS area from ISauIde^t alJ i|2015l) . adopting our same criteria on mass and size. Filled boxes plot the results using 
Sersic profiles, while open boxes are for the de Vaucouleurs profile (note that the results for the two profiles in the lowest redshift 
bin are superimposed). Right. Black open squares, solid lines and error bars plot KiDS number densities for candidate relic compacts, 
defined to have photometric formation redshift Zf > 2. Gray open squares and lines take into account possibl e systematics in the red shift 
determination. Long-dashed and solid green lines are for st ellar-like objec t s from BQ SS-DRIO with Zf > 2 l|Damianov et a n iioia and 
Mdyn > 6 X 10^*^ and 8 x 10^*^, respectively. The results from iTruiillo et alj ll2009l) and iTavlor et al ] <[20^) are shown he re as upper limits 
(see red and cyan points with arrows). The blue triangle is for the lower limit at 2 0 provided bv lTruiillo et all l|2014l). Orange and red 

lines plot abundances of relic galaxies from semi-analytical models (SAMs), based on Millennium N-body simulations llQuilis fc: Truiillq 
[ioia . Relics have been defined as system s whose ste l lar ma ss ha s increased since 2 = 2 to the present by less than 10% and 30%, 
respectively. Dashed and solid lines are for iGuo et adj[20]djj|}_andjGu o et alj (l2013h SAMs, respectively. Purple stars (connected by a 
line) are predictions from the lllustris simulations dWeilonT^ al.ll2015a^ . 


large-scale structure (i.e. the cosmic variance): they are cal- 
culated with the online CosmicVarianceCalculatoi0 tool 
llTrenti fc Stiavenill2008h . The cosmic variance increases the 
Poissonian error budget by ~ 5 — 30%. The total rela¬ 
tive error on abundances (i.e. number densities) amounts 
to ~ 25 - 30%. 

In Fig. |6]we plot the redshift evolution of the abundance 
of compact galaxies (left panel) and that for the subsample 
of systems with old photometric ages (i.e. formation red¬ 
shift Zf ^ 2; right panel). We consider these potentially old 
MSCGs, as candidate remnants of the compact ETGs found 
at high redshift by several studies (see Sec.[T]). We also re¬ 
determine the abundances by accounting for possible sys¬ 
tematics in the photometric redshifts (see Sec. 13.21 see grey 
symbols in Fig. EJ. We remind the reader that our sample 
starts to be incomplete at 2 ^0.5 (see vertical dashed lines). 

In the left panel of Fig. (6] we plot number densities 
for the whole sample of KiDS compacts, independent of the 
galaxy formation redshift. Our number densities are com¬ 
pared with estimates from IPamianov et al.l ll20l4l . for the 
number density of stellar-li ke objects havin g spectroscopic 
redshifts from BOSS-DRIO llAhn et al.ll20l4l and with three 
different cuts in total dynamical mass, Mdyn (> 5 x 10^°, 


^ http://casa.colorado.edu/~trenti/CosmicVariance.html 


8 X 10^°, and 1 0^^ Mp), respecti v ely; s ee Fig. El . Notice that 
objects in the lPamianov et al.l ll20l4l sample are not clas¬ 
sified according to either morphology or galaxy age, nor do 
they have an accurate estimate of the intrinsic R^. Hence, 
their selection may be missing compact systems that are ac¬ 
tually s£atially_£esolwdinMOSS-DR10. On the other hand, 
since iDamianov et al.l ||2014|) also include compacts with 
Re > 1.5 kpc and are selected with respect to Mdyn, rather 
than M*, it is not surprising, perhaps, that those abundance 
estimates are larger than ours. In fact, if we consider the 
abundance estimate of massive BO SS targets with R e < 1.5 
kpc and Mdyn > 8 x 10^° Mq from IDamianov et al.l ||2014| ) 
(see purple triangle in the left panel of Fig. El, the number 
density of compacts in the lowest redshift bin (2 ~ 0.25) is 
(1.7±0.5) X 10“® Mpc“®, in much better agreement with our 
density estimate of (1.6±0.4) x 10“® Mpc“® for the same red¬ 
shift bin. On the other hand, selecting galaxies with Re < 2.5 
kpc in KiDS, would lead to abundances ~ 0.8 dex larger than 
for Re < 1.5 kpc, still i n good agreement with estimates for 
stellar-like objects from lDamjanov et al.l (|2014| ). The yellow 
region in the left panel of Fig. El plots number densities f or 
galaxies in the COSMOS survey llDamianov et al]|2015al 'FI. 

® These data are kindly computed for us by 1. Damjanov (pri¬ 
vate comunication) by applying the same size and mass selection 
criteria as in the present work. 
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Remarkably, no evolution with redshift is found, for both 
KiDS and COSMOS samples, although, surprisingly, abun¬ 
dances for COSMOS (on average ~ 10~® Mpc~^) are about 
one o rder of magnitu d e larg er than our KiDS estimates. 
Since IPamianov et al.l (l2014l ) claim to find consistent den¬ 
sity estimates between COSMOS and BOSS (the latter hav¬ 
ing an area 4000 times larger than COSMOS), cosmic vari¬ 
ance seems not to be responsible for the above discrepancy. 
However, we cannot exclude that galaxy environment, which 
might be the actual dri ver of the number density of compact 
relics at 2: 0 (see, e.g. Tru jiUo et al.ll20l3 : IPoggianti et al.l 

l2013al : IValentinuzzi et al.ll20ldl . may be different for galax¬ 
ies in the KiDS-DR2 and COSMOS areas - an issue that 
will be addressed in forthcoming extensions of the present 
work. 

The results for photometrically old MSCGs (right panel of 
Fig. 16]) are fir st compared with estima tes for Zf > 2 compact 
galaxies from lPamianov et al.l ||2014| '). who selected samples 
with two cuts in Mdyn (> 6 x 10^° and 8 x 10^'’, respec¬ 
tively; see Fig. El). As for the left panel, there are differences 
in the selection criteria that make a direct comparison not 
trivial. However, the trends for both BOSS and KiDS sam¬ 
ples look qualitatively similar, with a very mild evolution 
with redshift in the range where our KiDS sample is com¬ 
plete. Moreover, a s noticed above, f o r BO SS compacts with 
Re < l.Skpc, the IPamianov et al.l ll2014h abundance esti¬ 
mates are fairly consistent with ours, within the uncertain¬ 
ties. 


Finally, we compare MSCGs number densities with 
predictions from sem i-analytical models 0 (SAMs). 
lOuilis fc Truiillol ll2013ll have determined the evolution of 
the abundance of compact galax ies from SAMs b a sed o n 
Millennium N-body simulations (IGuo et al.l 12011^ . 120131) . 
where relic compacts are defined as galaxies which have 
barely increased their stellar mass between 2: ~ 2 and 
2 ~ 0. Operatively, they selected from the merger tree 
those objects that have increased their mass since 2 = 2 
by less than 10% and 30% , respectively, i.e. galaxies 
whose mass at 2 ~ 2 is larger than 90% and 70% of the 
mass limit applied to select compacts. Our resu l ts are 
consistent in the lowest redshift bin with lGuo et al.l ll2013l l 
for simulated galaxies which have increased their mass at 
most by 10%, while we estimate lower densities than in 
simulations at higher redshift. However, similarly to what 
discussed for the comparison to BOSS estimates, theoretical 
predictions should be a ctually considered as upper limits, as 
lOuilis fc Truiillol ll2013ll did not apply any precise selection 
in size, since the resolution in the simulations does not 
allow reliable estimates of galaxy effective radii to be 
obtained. We also compare our findings to results from 
the hy drodynamical Illustris simu l ation of 
( 2015d ) (see also IWellons et al.l l2015bl ). 

( 2015al ) select 35 massive and compact galaxies at 2 = 2 
and follow their evolution to 2 = 0. Only 1 out of these 
35 systems evolves into a galaxy that satisfies our mass 
and size criteria at 2 ^ 0.5. This corresponds to a number 
density of 8.28 x 10“^Mpc“®, which is consistent with 


Wellnns_et_^ 

Wellons et al. 


^ We caution the reader that stellar masses and sizes are mea¬ 
sured in a different way between simulations and observations, 
hampering a straightforward comparison of the two. 


the abundances of compact galaxies that have accreted 
less than 10% of their final mass from the Millenium 
simulations. As a further caveat here, we point out that in 
our selection we adopt the same mass cutoff value at all 
redshifts, while simulations perform the mass selection at 
2 = 0, implying that at high redshifts, simulated galaxies 
with masses smaller than the 2 = 0 cutoff value are actually 
included in the analysis. This is another reason why one 
may expect that predicted number densities for compacts 
at high redshift are actually higher than the observed ones. 


At redshifts 2 < 0.2, in the left panel of Fig. El we 
see a l ack of candid a tes. T his seems to contrast the re¬ 
sults of iTruiillo et al.l ll2009ll who found 29 secure MSCGs at 
2 < 0.2 fulfilling our sa me criteria, all of them ha ving young 
ages < 4 Gyr (see also iFerre-Mateu et ^ I 2 OI 2 II. However, 
one sh ould notice that out of the 29 MSCGs of Trujillo et al.l 
ll2009li . only one is at redshift <0.1, still pointing to the ex¬ 
treme paucity of such systems in the nearby Univers e, con¬ 
sistent with our result. Similarly. iTavlor et al.l ll2010l) found 
one possible old MSCG at low redshift, using a more relaxed 
criterion for the size, than the one we adopt here. 


ISaulder et al.l ll2015l ) have found a sample of 76 com¬ 
pact galaxies from SDSS at 0.05 < 2 < 0.2, which resemble 
quiescent galaxies at high-z, i.e. systems with small effective 
radii and large velocity dispersions. If we consider their com¬ 
pacts with Re < 1.5 kpc and M* > 8 x 10^'^ Mq, 1 galaxy 
at 2 < 0.1 and 6 galaxies at 2 > 0.1 are left when the effec¬ 
tive radius from a de Vaucouleurs profile is used. Instead, 
these numbers change to 1 galaxy at 2 < 0.1 and only 1 at 
2 > 0.1 if a Sersic profile is fitted. These numbers correspond 
to abundances of 2.4 x 10“® Mpc“® in the redshift range 
0.05 < 2 < 0.1, and 2 x 10~® Mpc“® and 3.3 x 10~® Mpc“^ 
at 0.1 < 2 < 0.2, if de Vaucouleurs or Sersic profile are 
fitted, respectively. As mentioned in Sec. [T] these Hndings 
seem to trouble the current hierarchical paradigm of galaxy 
formation, where some relic syste ms at 2 ~ 0 are actually 
expected to be found. In contrast, IPoggianti et al.l (l2013ah 
have found 4 galaxies fulfilling our same criteria (correspond¬ 
ing to 1.4% of their sample galaxies with masses larger than 
8 X 10^'^ Mq), and all of these galaxies are old, with mass- 
weighted ages older than 8 Gyr. These numbers translate 
into a very l a rge a bundance of ~ 10~®Mpc“®. Recently, 
iTruiillo et al.l (l2014h added a new brick to the story, finding 
one relic compact in the Perseus cluster, i.e. NGC 1277, 
reconciling the observations at 2 ^ 0.2 with predictions 
from simulations. All these results point to an overabun¬ 
dance of MSCGs in dense cluster regions, that are actually 
under-represented over the area currently mapped by KiDS. 


As discussed in Sec. E21 the sample of MSCGs would be 
reduced significantly in size by systematics in the redshift 
determination. Fig. El“corrects” the abundance estimates for 
the possible systematics in redshift. The corrected abun¬ 
dances (see gray boxes and lines in Fig. El should be seen 
as our current lower limit on number densities of SMCGs 
in KiDS. This issue will be addressed by on-going spectro¬ 
scopic follow-up of MSCGs, and discussed in details in a future 
paper. 
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4 CONCLUSIONS 


Thanks to the large area covered, high image quality, excel¬ 
lent spatial resolution and seeing, the Kilo Degree Survey 
(KIDS) provides a unique dataset to study the properties of 
super-compact massive galaxies (MSCGs) - a family of sys¬ 
tems which plays a key role into our understanding of galaxy 
formation and evolution. 

In this paper, we present a sample of candidate MSCGs , 
based on 156 sq. deg. of KiDS, in four optical bands (w, g, r 
and i). We start from a sample of ~ 0.4 million galaxies 
with high-S/N, measured photometry and structural pa¬ 
rameters in all the four bands. For a subsample of these 
galaxies, we have used the KiDS photometry to estimate: 1) 
photometric redshif t s based on machine learning techniques 
jBrescia et al.l [20141 : ICavuoti et al.ll2015l ): 2) structural pa¬ 
rame ters using the software 2DPHOT llLa Barbera et all 
l2008tl : 3) stellar masses, fitting colours with SPS model pre¬ 
dictions. The resulting sample is > 90% complete down to 
an r-band magnitude ~ 21, and down to a stellar mass of 
3 — 5 X 10^° Mq , up to a redshift 2 ~ 0.5. We select the most 
massive (M* > 8 x 10^° ) and most compact {Re < 1.5 kpc) 
galaxies with (photometric) redshift 2 < 0.7. We remove 
star contaminants by performing a visual inspection of the 
final sample of candidates and then, for galaxies with avail¬ 
able near-IR photometry from VIKING-DRl, we combine 
optical-fNIR photometry to reduce the fraction of contami¬ 
nants. 


The final sample consists of 92 compact candidates, 
with a number density of ~ 0.9 compact galaxies per square 
degree, at 2 < 0.7. Nine candidates have spectroscopic infor¬ 
mation from SDSS and GAMA surveys, that is used to assess 
the systematics in the redshift determination of our sample. 
On average, compact galaxies have negative colour gradients 
which are simila r to the ones for normal passively evolv¬ 
ing galaxies (e.g., Tamura et al.ll2005 : Tamura fc Ohtall2000l : 


Tc^ora et al.l l20ld : iTortora et al.l l2013l : iLa Barbera et ahl 
2OI2II . However, the variety of gradients, to be confirmed 


with a spectroscopic follow-up of the present sample, seems 
to suggest that comp act galaxies formed u nder a wide range 
of initial conditions (iGargiulo et al.l[2012l i. We also discuss 
the evolution with redshift of the number density of com¬ 
pact systems, and in particular that for the oldest galaxies, 
which are possibly remnants of the high-redshift (2/ ^2) 
compact population detected by many studies. Remarkably, 
we do not find any MSCGs candidate at 2 < 0.2. This find¬ 
ing, which is consistent with the recent studies, might be 
related to the effect of galaxy environment on the abun¬ 
dance of compact systems (see Sec. 1). Although obser¬ 
vational studies at intermediate- and high-z do not point 
to a clear picture of how e nvironment affects galaxy sizes 
(see iDamianov et al.ll2015bl and references therein), recent 
cosmological simulations predict a larger fraction of mas- 
sive compact system s in high- than in low-density regions 
llStringer et al.|[2015f l. This prediction is supported by the 
fact that NGC1277 - the only well-characterized compact, 
massive, ETG at 2 ~ 0 - res ides in the core reg ion of a 
nearby rich cluster of galaxies ^ Trujillo et al.ll20l4: see also 
IValenitinuzzi et ^l201(ll : IPoggianti et al. 2013a 1. Thus, the 
absence of compact galaxies at 2 < 0.2 could be related to 
the smaller fraction of dense structures in the area currently 
mapped by KiDS. 


We plan to have 10 times more compact candidates at 
the end of the KiDS survey, when all the 1500 sq. deg. will 
have been observed. Only nine of our candidates have a 
spectroscopic coverage at the moment. Thus, a spectroscopic 
follow-up is actually necessary, to fully validate and charac¬ 
terize our sample, and determine accurate number densities 
to be compared with theoretical expectations. 
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